Gamma-ray Constraints on Dark Matter Annihilation into Charged Particles 
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Dark matter annihilation into charged particles is necessarily accompanied by gamma rays, pro- 
duced via radiative corrections. Internal bremsstrahlung from the final state particles can produce 
hard gamma rays up to the dark matter mass, with an approximately model-independent spectrum. 
Focusing on annihilation into electrons, we compute robust upper bounds on the dark matter self 
annihilation cross section ((7A«)e+e- using gamma-ray data from the Milky Way spanning a wide 



range of energies, 



10" 



10 GeV. We also compute corresponding bounds for the other charged 



leptons. We make conservative assumptions about the astrophysical inputs, and demonstrate how 
our derived bounds would be strengthened if stronger assumptions about these inputs are adopted. 
The fraction of hard gamma rays near the end point accompanying annihilation to e'^e~ is only 
a factor of < lO'^ lower than for annihilation directly to monoenergetic gamma rays. The bound 
on {o-Au)e+e- Is thus Weaker than that for {ctav)~,~i by this same factor. The upper bounds on 
the annihilation cross sections to charged leptons are compared with an upper bound on the total 
annihilation cross section defined by neutrinos. 

PACS numbers: 95.35.-(-d, 95.85.Pw, 98.70.Vc, 98.62.Gq 
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I. INTRODUCTION 

A wealth of observational evidence attests to the exis- 
tence of dark matter (DM) (see, e.g., Refs. [^-Q fo'" 
views). However, despite the fact that DM accounts for a 
large fraction of the total energy density of the Universe, 
it has evaded direct detection and its particle properties 
remain unknown. 

If dark matter is a thermal relic, it should have a small, 
but non-negligible, self annihilation cross section. This 
allows a means of indirect detection of DM, as annihila- 
tions in the Universe today may generate an observable 
flux of products. For thermal relic DM, a total anni- 
hilation cross section of (ctav) ~ 3 x 10~^®cm'^ s~^ is 
needed in order to obtain the observed relic abundance 
of r^DM — 0.3. (We shall work with the quantity {(Tav), 
which is the product of the annihilation cross section and 
relative velocity, averaged over the dark matter velocity 
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-3, 



If DM is 



distribution. In the Milky Way, Vrr. 

not a thermal relic, e.g., Refs. [3-|6|, even larger annihi- 
lation cross sections are possible. There have also been a 
number of recent proposals in which {uav) is enhanced 
at low velocity |7h9|. Il4l - [2(]j (i.e., in galactic halos) while 
still satisfying the thermal relic constraints. (Note that 
these scenarios are subject to constraints arising from 
annihilations in the first collapsed structures [HI.) 

While the total cross section may be related to the 
relic abundance, the branching ratios to any particular 
final states are model-dependent. If we assume the DM 
is the lightest new particle beyond those in the stan- 
dard model (SM), then the branching ratio to SM final 
states must be 100%. Annihilations must then produce 
fiuxes of detectable particles emanating from regions of 
DM concentration, with signals such as gamma rays, mi- 
crowaves, neutrinos, and positrons being of particular in- 
terest. There exist general constraints on the total anni- 
hilation cross-section, based upon unitarity [2^ . [23| and 



the requirement that annihilations not significantly alter 
halo density profiles [1]. A strong bound is placed on 
the total annihilation cross section by assuming that the 
branching ratio to neutrinos, the least detectable final 
state, is 100% [H,!!^. 

We focus here on the process XX ~^ e+e" in which 
DM annihilates to an electron-positron pair, though we 
shall also report results for the other charged leptons. Al- 
though the branching ratio to this particular final state 
is model-dependent, it is a significant channel in a wide 
range of models. For example, while annihilation to 
fermions is helicity suppressed in supersymmetric mod- 
els, Kaluza-Klein DM features large (unsuppressed) an- 
nihilation rates to leptons p^| - [l^ , as does the Dirac DM 
model of Ref . [l3| . Numerous authors have recently pro- 
posed models in which annihilation to charged leptons 
is significantly enhanced [l4'J20'| , making upper limits on 
the cross section to these annihilation products partic- 
ularly interesting. In addition, various other SM final 
states, such as W~^W~ and ZZ, produce via their 
decays and hence a flux of charged leptons is of generic 
interest in a large variety of DM models. 

Several techniques may be used to constrain the pro- 
duction of e~^e~' within galactic halos, all of which rely 
on the fact the charged particles inevitably produce pho- 
tons. Signals considered include gamma rays, x-rays, mi- 
crowaves and radio waves. Photons are produced by the 
various energy loss processes that charged particles un- 
dergo in a galactic halo, examples of which include syn- 
chrotron radiation due to the propagation of in galac- 
tic magnetic fields, and inverse Compton scattering of 
electrons from interstellar radiation fields, e.g., |26l - l3a |. 
The drawback of these techniques is a significant depen- 
dence on astrophysical inputs, some of which are poorly 
known. Uncertainties in magnetic field strengths, radia- 
tion backgrounds, and electron diffusion scales all enter 
the calculations in an involved fashion. 
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Charged particles also produce photons via electro- 
magnetic radiative corrections |36l - l44| . The lowest order 
dark matter annihilation process XX ~^ e^e^ is necessar- 
ily accompanied by the radiative correction xX ~^ e~^e~"f. 
This is an internal bremsstrahlung (IB) process, mean- 
ing that the photon arises at the Feynman diagram level 
and is not due to interaction of charged particles in a 
medium. Importantly, for a given annihilation cross sec- 
tion, {aAv)e+e-i the accompanying flux of IB photons 
can be determined without knowledge of the new under- 
lying particle physics which mediates the DM annihila- 
tion. Moreover, IB suffers none the drawbacks of the 
competing inverse Compton and synchrotron techniques 
outlined above. While inverse Compton and synchrotron 
fluxes are dependent on conditions of the astrophysical 
environment, the IB flux is always present and its nor- 
malization and spectrum are predetermined. 

In this paper we use IB emission to derive robust up- 
per limits on the dark matter annihilation cross section to 
electron-positron pairs {(JAv)e+e- over a wide DM mass 
range spanning ^ 10~^ — 10'' GeV. We calculate DM 
annihilation fluxes produced in the galactic halo, and 
compare with the gamma-ray backgrounds reported by 
COMPTEL, EGRET and H.E.S.S. We also look at data 
for the M31 (Andromeda) galaxy, to fill a gap between 
the energy ranges covered by EGRET and H.E.S.S. We 
explicitly demonstrate how our limits vary according to 
the assumed DM halo profile (our one source of uncer- 
tainty) and also compare our limit on the annihilation 
cross section to e"'"e^ with corresponding bounds on the 
77 and vu final states. 




FIG. 

XX - 
GeV. 



E [GeV] 



1: Internal bremsstrahlung gamma-ray spectra per 
• e+e" annihilation, for = 100, 200, 500, and 1000 



This spectrum is shown in Fig. [T] for various choices of 
the DM mass, where a sharp edge in the spectrum at 
E — nix evident. 

Note that we consider only radiation from the fi- 
nal state particles, and not from any internal propa- 
gators. In some supersymmetric scenarios in which s- 
channel annihilation to fermions is helicity suppressed, 
bremsstrahlung from internal propagators can be par- 
ticularly important as it can circumvent this suppres- 
sion pt]| - |43| . We do not consider these model-dependent 
processes. Note, however, that the presence of such emis- 
sion would only increase the gamma-ray flux we calculate, 
and hence strengthen the cross section limits derived. 



II. INTERNAL BREMSSTRAHLUNG 

If DM annihilates to produce charged particles, the 
lowest order processes will always be subject to elec- 
tromagnetic radiative corrections, resulting in the pro- 
duction of real photons. In particular, the annihila- 
tion XX ~^ e~^e~ will be accompanied by the internal 
bremsstrahlung process xX ~^ e~^e~j. A photon may be 
emitted from either the final state e"*" or e~ , with a cross- 



section proportional to a ~ 1/137. See Refs. [36l - l38l . l45l| 
for a detailed discussion. To a good approximation, the 
differential cross section for xX ^ e+e^7 is 
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where E is the photon energy, s = 4m^, s' = 4TO^(m^ — 
E), and Ctot is the tree-level cross section for xX e+e^. 
Note that crtot factors from the IB cross-section. This 
important feature implies that the IB spectrum is in- 
dependent of the unknown physics which mediates the 
lowest order annihilation process. The photon spectrum 
por XX ~^ e"^e~ annihilation is therefore given by 
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III. ANNIHILATION IN DARK MATTER 
HALOS 

The rate at which dark matter annihilates is propor- 
tional to the square of the dark matter number density, 
= p/rrix- However, there are considerable uncertain- 
ties in the dark matter density profile of the galactic halo, 
particularly in the central region where the density, and 
hence the annihilation rate, is largest. To deal with these 
uncertainties, we make conservative assumptions about 
the DM density profile, and show how our results vary if 
less conservative assumptions were used. 

A standard parameterization of a dark matter halo 
density profile is 



p{r) = 



{r/rsV [1 + [r/rs) 



a-\(l3-l)/a 



(3) 



The Kravtsov [401, Navarro- Frenk- White (NEW) ^ and 
Moore [48] profiles are defined by the values (a,/3,7) 
shown in Table ID The Milky Way values for rs, the 
scale radius, and p(i?sc), the DM matter density at the 
solar circle i?sc = 8.5 kpc, are also given in Table H] The 
normalization po is then fixed by p(i?sc) and r^. For large 
radii f > r^, all three profiles scale with radius as 
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TABLE I: Values of the parameters {a, P, 7) which define the 
Kravtsov, NFW and Moore halo profiles. The values of the 
scale radius rs in [kpc] , and halo density normalization at the 
solar circle p(-Rsc) in [GeV cm~^], are specific to the Milky 
Way 



where dNj/dE is the gamma-ray spectrum per annihi- 
lation. For the IB emission associated with annihila- 
tion to e^, we must replace {uav) with {(jA'v)e+e- = 
{gav) X Br(e+e^), while dN^/dE is given by Eq. [51 
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and are normalized such that they coincide closely. How- 
ever, the profiles diverge for small radii, scaling as l/r°-^, 
l/r and for the Kravtsov, NFW and Moore pro- 

files respectively. The steep Moore profile thus features a 
greatly enhanced density near the Galactic center, com- 
pared to the relatively flat Kravtsov profile; the NFW 
profile falls between the two. 

Uncertainties in the density of the DM halo translate 
into uncertainties in the DM annihilation rate. A detailed 
discussion of the dependence of DM annihilation signals 
on the choice of density profile is given in Ref. [25]. While 
the uncertainties are mild for large angular regions of the 
galaxy, they scale several orders of magnitude for small 
angular scales close to the Galactic center. In order to 
place conservative upper limits on the DM annihilation 
cross section, we focus on the profile with the smallest 
dark matter density, namely the Kravtsov profile. For 
the NFW and Moore profiles, smaller values of {<Jav) are 
needed to reproduce the same flux, and hence lead to 
stronger (less conservative) limits. 

We now calculate the gamma-ray fiux for annihilations 
in the galactic halo. For an observation direction at angle 
-0 with respect to the Galactic center, the integral of the 
square of the dark matter density along the line of sight 
is given by 



JW = Jo 



^ """^ {^^Rl-^^Rsc cos^ + e^'^ di , (4) 



where Jq = l/[8.5kpcx (0.3 GeV cm"^)^] is an arbitrary 
normalization constant used to make J{ip) dimension- 
less, and which cancels in our final expression for the 
gamma-ray flux. We then define the average of ^7(0) 
over an observation region of solid angle Ail as 



2n /"^ 

Jau = J J {ip) sin ipdil; 



(5) 



Values of J'(V') and Jau are given -0 in Fig. 2 of Ref. |25j . 
With these definitions, the gamma-ray flux per steradian 
due to DM annihilation in an observation region of an- 
gular size Af2 is 



rf^T- _ {(JAv) Jaq. dN^ 
dE ~ 2 ATTmpo dE 



(6) 



IV. ANALYSIS OF ANNIHILATION FLUX 
A. Analysis technique 

Our analysis technique is similar to that followed in 
Mack et al. Ea- We use galactic gamma-ray data from 
COMPTEL 5^, EGRET |5ll and H.E.S.S. [sl, together 
spanning the broad energy range 10~^ — 10'* GeV. As 
there is a small gap between the energy ranges covered 
by EGRET and H.E.S.S., we use the observations of the 
M31 (Andromeda) galaxy made by CELESTE ^] to 
calculate constraints for this energy interval. 

The galactic gamma-ray background measurements re- 
ported by COMPTEL, EGRET and H.E.S.S. are given 
in approximately log-spaced energy intervals, with energy 
bins of size ranging from AlogE ~ 0.2 — 0.6. We cal- 
culate the IB gamma-ray flux for the observation regions 
viewed by these experiments, using the methods outlined 
above, and compare with the observational data. Upper 
limits on {aAv)e+e- are determined by requiring that the 
IB flux due to DM annihilation be lower than 100% of 
the observed gamma-ray background flux in each of the 
experimental energy bins. Given that a large fraction of 
the observed gamma-ray background is likely to be as- 
trophysical in origin, and not due to DM annihilation, 
taking the total background flux to be an upper limit on 
the DM annihilation signal is an extremely conservative 
approach. 

For each energy bin, we take the DM mass to be equal 
to the upper energy limit of the bin, and integrate the 
IB flux over the width of the bin. Figure [5] shows the 
quantity 
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which is the number of photons per annihilation as a 
function of bin size, for IB emission from e*, and m-^ = 
1000 GeV. [Since enters Eq. [T] via a logarithm, vari- 
ation with is only very mild.] This indicates how 
the size of the energy bins affects our results. For com- 
parison, the number of photons per annihilation for the 
process xX ~^ 11 (fo^ which dN^/dE — 2S{m^ " E)) 
is also shown. Note that the IB cross section is propor- 
tional to 1/E, so for sufficiently low photon energy the 
IB probability becomes large and one must account for 
multiple photon emission. However, we are not working 
in this regime, and in fact obtain our limits using the fiux 
near the endpoint of the spectrum. Despite the fact that 
the IB flux is small, we have enough hard gamma rays 
near the endpoint to result in strong bounds. For typical 
parameters, the IB flux per annihilation to is smaller 



4 



10 E 



Typical 
Bin Size 




0.001 



FIG. 2: Number of gamma rays per DM annihilation 
{f^^_ dE dNj/dE) as a function of the lower limit of in- 
tegration, for IB emission from XX ~^ e'^e" (solid Hne). A 
typical bin size used in the analysis is shown. The DM mass 
used is 1000 GeV; variation with m-f^ is very small. Shown 
for comparison is the number of photons per annihilation for 
the process xX ~^ 77 (dashed line) in which the photons are 
always at the endpoint. 



than the photon flux per annihilation for xX — 77 by a 
factor of 10^. This is expected, given that the IB cross 
section is suppressed by a factor of a with respect to the 
tree-level DM annihilation process. 



Observational data 
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FIG. 3: Upper limit on {gv)^+^- as a function of DM mass 
for the Kravtsov (solid line), NFW (dashed line) and Moore 
(dotted-dashed line) profiles. 



a 2-(T upper limit on the flux of gamma rays from M31 
between 50 and 700 GeV of around 10"^*^ photons cm^^ 
s~^ was reported 'sT] . We compare this flux with the IB 



signal calculated for an energy bin of size 10" 
This is extremely conservative, as we are constraining the 
cross section by requiring that the annihilation flux in a 
small bin be less than or equal to the observed flux in 
a much larger bin. The DM density profile of M31 is 
less well constrained than that of the Milky Way. As 
in Ref. [4^, we model the M31 halo using the Kravtsov 
profile for the Milky Way, which yields Jao x Afi ~ 
2 X IQ-^. 



COMPTEL and EGRET are telescopes aboard the 
Compton Gamma-Ray Observatory. Refs. [13, HE] 
present flux data from these telescopes over an energy 
range of 1 MeV to 100 GeV between them, with a ob- 
servation region of —30° < I < 30° in Galactic longi- 
tude, and —5° < 5 < 5° in Galactic latitude. We cal- 
culate Jao as if the DM annihilation signal were from 
a circular region of -0 = 30°. This is a conservative ap- 
proach, as the average annihilation flux per steradian for 
the circular region is lower than for the rectangular re- 
gion which was actually observed, and yields values of 
Jau = (13,28, 100) for the Kravtsov, NFW and Moore 
profiles respectively. 

The H.E.S.S. observations reported in Ref. [s^ cover 
the relatively small angular region —0.8° < I < 0.8°, 
—0.3° < 6 < 0.3°, in Galactic coordinates, over an en- 
ergy range of 300 to 15 000 GeV. The data include a 
background subtraction which we must take into account 
when calculating Jaq] see Mack et al. [49| for details. 
This leads to Jaq = (3,850,50 000) for the Kravtsov, 
NFW and Moore profiles respectively. 

CELESTE viewed the M31 galaxy with an observation 
region of angular radius 0.29°. No signal was seen, and 



V. DISCUSSION 

In Fig. [3]we show the upper limits on {aA'v)e+e- as a 
function of DM mass, using the observational data de- 
scribed above. We give the Galactic genter results for 
the conservative Kravtsov profile, the more commonly 
adopted NFW profile, and the steep Moore profile. For 
the CELESTE observation of M31, differences between 
these profiles are expected to have a modest effect on the 
results, as a large portion of the galaxy is within the field 
of view; in Fig. [3] we show the CELESTE constraint using 
only the Kravtsov profile. As previously discussed, while 
the Kravtsov, NFW and Moore profiles diverge towards 
the center of the Galaxy, they are similar at large radii. 
As the EGRET and COMPTEL observations encompass 
relatively large angular scales, the density profile changes 
have a modest effect. On the other hand, the H.E.S.S. 
constraints correspond to a much smaller angular region 
toward the Galactic center, and vary by orders of mag- 
nitude depending on the profile adopted. (See Ref. [23 
for a full discussion of the differences between the pro- 
files for different angular regions.) To be conservative, 
we do not consider the possibility that DM annihilation 
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FIG. 4: Upper limits on the partial cross sections Br{ii) x 
{o"v)totai for various final states ii — e^e" (solid black line; la- 
beled), fj,^ (thick dashed line; labeled), t'^t~ (thick dashed 
line; labeled), 77 (red line; labeled), and j/f (blue line; la- 
beled) , using the conservative Kravtsov profile. Each of these 
partial cross-section limits is independent, with no relation- 
ship assumed between the branching ratios to particular final 
states. Also shown are the KKT (thin dashed line) and unitar- 
ity (thin dotted-dashed line) limits on the total cross section 
described in the text, and the cross section for thermal relic 
DM (natural scale). The 77 and 1^1/ limits are taken from 
Ref. gi] and Ref. gl], respectively. 



rates are enhanced due to substructure in the halo, e.g.. 



, or mini-spikes around intermediate-mass 
6^ IgsI : such enhanced annihilation signals 
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black holes 

would result in stronger upper bounds on the cross sec- 
tion. 

We can estimate the way that the cross-section bounds 
scale with the DM mass. The spectrum of the galactic 
gamma-ray background falls off with energy as (i$ /dE ^ 
E^" where, e.g., a is shghtly larger than 2 in the EGRET 
and HESS energy ranges. The IB signal scales approx- 
imately as d^/dE ~ (aAv)m~^E~^ [where accounting 
the full energy dependence in Eq. [T] has only a small ef- 
fect on this scaling]. Given this scaling with E, it is clear 
that the strongest constraints arise from the endpoint of 
the spectrum. In addition, small bin size is optimal for 
obtaining strong constraints. If we integrate the flux over 
an energy bin of width xm^ to , the IB flux within the 
energy bin is proportional to {crAv)'m^^ (ignoring loga- 
rithmic corrections) while the background flux is propor- 

The cross-section limits 
and thus rise as 



tional to m^"+^ (for fixed x) 

(CTAV) 



X 

then scale with 
mv if a ~ 2. 



mv as 



m" 



In Fig.|3]we show the upper bounds on the annihilation 



cross sections into e+e~, iJ.^fJ,~, and r+r", based upon 
the IB emission from each final state (all use the conserva- 
tive Kravtsov profile) . As the mass of the charged lepton 
increases, the rate of IB emission decreases and thus the 
upper bounds on the cross sections become weaker. How- 
ever, as the IB spectrum depends only logarithmically on 
the charged lepton mass, this effect is mild, particularly 
for large DM mass. We have used Eq. [l]to calculate the 
IB flux over the entire mass range, and present limits 
which range from high down to just above threshold. 
We expect modifications to Eq. [T] in the limit that the 
charged leptons are nonrelativistic, but this will only af- 
fect a small mass range close to threshold. In the case of 
annihilation to , we have not considered the gamma 
rays that arise from hadronic decay modes of the r lep- 
tons (see, e.g., Ref. [38|) which form a broad spectrum 
centered on the pion mass. 

Using Fig. 21 we may compare the limit on the anni- 
hilation cross section into e~^e~ with that for 77, taken 
from Ref. [1^. As anticipated, the bound on (cru)g+g- 
from IB is weaker than the bound on {av)-y-y by a fac- 
tor of ~ 10~^ ~ a. This difference can be understood 
by comparing the number of hard gammas near the end- 
point, E^ = m^. For annihilation to 77, there are always 
two monoenergetic gamma rays at the endpoint. The IB 
spectrum, integrated over a typical bin width (see Fig. [5]) 
results in a flux that is smaller than this by only a factor 
of less than 100. That IB provides a limit this close to 
the ideal 77 channel illustrates the importance of the IB 
technique. 

We stress that our IB limits apply to the partial DM 
annihilation cross sections, {(tav)i+i- , rather than the to- 
tal annihilation cross section, (cfav)- The two are related 
via {(JAv)i+i- — Br{l^l^){aAv). While the branching 
ratios are entirely dependent on the choice of DM model, 
there are many scenarios which feature large branching 
ratios for direct annihilation to leptons. such as Kaluza- 
Klein DM, or the recent models of Refs. |14l - [20l |. Our con- 
straints can be readily applied to any particular model, 
simply by dividing the (aAv)i+i- bounds by the relevant 
branching ratios. 

Also shown in Fig. 2] are a number of upper bounds 
on the total annihilation cross-section, based upon uni- 
tarity [22|, [2^ and the requirement that DM annihilation 
not significantly alter halo density profiles [Kaplinghat- 
Knox- Turner (KKT) Q]. Note that the limit on anni- 
hilation to Di/ also defines a strong bound on the total 
DM annihilation cross section [23| . If we assume annihi- 
lation to only standard model final states, a conservative 
bound on the total cross section is obtained by assum- 
ing the branching ratio to neutrinos (the least detectable 
final state) is 100%. Any other assumption would lead 
to appreciable fluxes of gamma rays and hence be more 
strongly constrained. Dark matter annihilation into neu- 
trinos was examined in Refs. [3l and [2^. The upper 
bound on ((TA'y)pi/ shown in Fig. U] (taken from (2^) shows 
that the neutrino constraints are quite strong, particu- 
larly for large DM mass. (Reference [63| extended the 



6 



DM annihilation limits to lower masses, while Ref. [65| 
included substructure enhancement. Reference [6l| de- 
rived analogous limits on the DM decay rate.) In fact, 
due to electroweak bremsstrahlung (radiation of W and 
Z bosons, rather than photons) neutrinos are also in- 
evitably accompanied by photons [67l - [69j |. However, this 
results in bounds that are comparable to, or weaker than, 
those obtained directly from neutrinos. 

If DM annihilates to e^, photons will be produced not 
only by IB, but also by energy loss processes including 
inverse Compton scattering and synchrotron radiation. 
In particular, radio wavelength signals produced via syn- 
chrotron emission have been the focus of much recent 
attention, e.g., (30l - [33 |. However, the intensity of syn- 
chrotron radiation depends on a number of uncertain as- 
trophysical parameters, such as magnetic field strength, 
radiation field intensities, and electron diffusion scales. 
By contrast, IB is free of these astrophysical uncertain- 
ties, and has a fixed spectrum and normalization. An- 
other key difference is the energy of the photons. Syn- 
chrotron radiation produces generally low energy pho- 
tons, while IB provides some hard gamma rays near the 
endpoint. Since the background flux falls off with energy, 
these hard gamma rays are extremely useful. The sharp 
edge in the IB spectrum at £' = TOjj, can be used to diag- 
nose the DM mass; this is not possible with synchrotron 
radiation. 

Nonetheless, it useful to take the synchrotron-based 
cross-section bounds as a reference point to compare 
with our IB-based bounds. Our conservative IB bound 
on {aAv)e+e- is comparable to conservative bounds 
on {aAv)w+w- obtained from synchrotron radiation. 
For example, Ref. [sJl obtains {aAv)w+w- ^ 4 x 
lO-^^cm^s^i (4 X 10-23cm3s-i) at 100 GeV (1 TeV) 
assuming an NFW profile and conservative magnetic 
field choices (lower panel of Fig. 6 in Ref. jS-i]). This 
is to be compared with our IB result of {aAv)e+e- ^ 
2 X 10-23cm3s-i (4 X IQ-^^cm^s-i) at 100 GeV (1 TeV), 
again assuming an NFW profile. (The results of Ref. [3J| 
are very similar to those of Ref. '32'] , though weaker than 
those of Ref. [31], in which less conservative assumptions 
were made.) Note that these synchrotron studies assume 
annihilation to W'^W~ (or qq) which then decay to elec- 
trons, rather than direct annihilation to e+e" . Therefore, 
these electrons are not at the DM mass, and have instead 
a broad distribution of energies centered on the W mass. 
The synchrotron analyses in [3l|, IH, HJ] thus serve only 
as an interesting reference point for our work, and not as 
a direct comparison. 

References (T^ - [20| have recently proposed models in 
which DM annihilates directly to charged leptons, with 
cross sections well above that expected for a thermal 
relic. This may account for anomalies in cosmic ray 



spectra from PAMELA, HEAT and ATIC, gamma-ray 
measurements from EGRET, and microwave signals from 
the Wilkinson Microwave Anisotropy Probe, all of which 
seem to require more electrons and positrons than can 
be explained otherwise. Our bounds on {aAv)i+i- will 
directly constrain the allowed parameter space for these 
types of DM models. 

We expect the sensitivity of the IB bounds to be im- 
proved by forthcoming data from the Fermi-GLAST ex- 
periment [70I - I72I . Improved point source subtraction en- 
abling the diffuse background to be reduced, together 
with better energy and angular resolution and high statis- 
tics measurements, will enable stronger limits to be 
placed on all DM annihilation processes that produce 
gamma rays in the measured energy range. 



VI. CONCLUSIONS 

Dark matter annihilation into charged particles will 
necessarily be accompanied by gamma rays. Internal 
bremsstrahlung from final state charged particles can 
produce hard gamma rays, close to the endpoint defined 
by — m-^, with an approximately model-independent 
spectrum. Using galactic gamma-ray data, we have cal- 
culated upper limits on the dark matter annihilation 
cross section to e+e^ and other charged leptons. We have 
made conservative assumptions about the astrophysical 
inputs, and demonstrated how our derived bounds would 
be strengthened if the galactic halo has a steeper density 
profile than assumed. The upper bound on the annihila- 
tion cross section into e^e^ is weaker than that for the 
ideal 77 final state by only a factor of < 10^. For a 
wide range of masses, our upper bound on {aAv)e+e- is 
stronger than the bound on the total cross section de- 
fined by neutrinos, the least detectable final state. Com- 
pared with recent constraints on DM annihilation cross 
sections based upon synchrotron radiation, the inter- 
nal bremsstrahlung constraints on {aAv)e+e- are broadly 
comparable in strength. However, synchroton emission 
depends strongly on poorly known astrophysical inputs, 
such as galactic magnetic field strengths. In comparison, 
the normalization and spectrum of IB radiation is fixed, 
independent of any astrophysical inputs, and is thus an 
extremely clean technique. 
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